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AN EXPERIMENTAL INVESTIGATION AND 
LIMIT ANALYSIS OF NET AREA IN TENSION* 


W. G. Brady** and D. C. Drucker*** (Brown University) 


Abstract 


The old question of the weakening effect of holes in tension members is 
re-examined from the viewpoint of plastic action. Using the maximum 
shearing stress yield criterion and the Mises criterion, upper and lower 
bounds for the limit load are computed for thin sheets with six configurations 
of staggered rows of holes. The limit load is defined as the force required 
for unrestricted flow across the entire member. Comparison is made with 
current design specifications, and with the results of a large number of 

tests performed on hot-rolled steel and on 61S-T6 aluminum alloy sheet. 
Plugs were inserted in the holes in a few of the tests. It is found that the 
limit loads, as determined by test, confirm the main features of the analysis, 
and lie somewhat above the ranges of load computed with the maximum 
shearing stress criterion. The s2/4g design rule is shown to correspond to 
an approximate upper bound, and as such, may be appreciably on the unsafe 
side for configurations of open or plugged holes, although probably reason- 
able for riveted tension joints. 


Introduction 


The recently developed theorems for limit analysis!2 offer an opportuni- 
ty to check current design specifications for net area in tension members, 
as well as to determine how closely the limit loads of selected configurations 
of rows of holes in steel and aluminum alloy sheets, as determined by test, 
agree with idealized theory. The limit load is here defined as the force 
required for unrestricted plastic flow to occur across the entire section and 
is therefore above the load at which plastic deformation first occurs. At 
this stage, either localized or overall plastic deformation is considerable so 
that the limit load may provide a reasonable value upon which to base the 
working load. 


*The results presented in this paper were obtained in the course of 
research sponsored by the Office of Naval Research under Contract 
N7onr 35801. 

**Research Assistant in Applied Mathematics. 

***Chairman, Division of Engineering. 

1. “The Safety Factor of an Elastic- Plastic Body in Plane Strain”, by 
D. C. Drucker, H. J. Greenberg, and W. Prager, Journal of Applied 
Mechanics, vol. 18, 1951, pp. 371-378. 

2. “Extended Limit Design Theorems for Continuous Media”, by D. C. 

Drucker, W. Prager, and H. J. Greenberg, Quarterly of Applied 

Mathematics, vol. 9, 1952, pp. 381-389. 
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Several design specifications have been proposed and used for determin- 
ing net area in tension members with staggered rows of holes. For one row 
of holes normal to the length of the member, subtracting the hole diameters 
from the gross width of the member is a rational rule for determining net 
width under the assumption of uniform stress across the section after plastic 
flow has occurred. For holes that are ataggered across the normal section, 
the implications of an assumption of a uniformity of stress due to plastic 
deformation are not so obvious. Consider two staggered holes as in Fig. 1. 
There are two extreme possibilities: one, to take the net area across each 
normal section, i.e. W-d, as net width, and the other, to take out both holes, 
i.e. W-2d. Here, W is the gross width of the piece, and d the diameter of 
each hole. Intuitively, it would seem that the first case would be unsafe, and 
the second too conservative. The correct net width probably lies between 
these extremes, and, moreover, it seems likely that the correct net width is 
a function of the stagger and gage, stagger being defined as the spacing longi- 
tudinally between the centers of adjacent holes, and gage the transverse 
spacing between centers of two adjacent holes. 

Many suggestions were advanced for design rules?»4 and in 1915 T. A. 
Smith5, on the basis of Cochrane’s earlier work, published the following 
formula: 


where 


x = the fractional part of a hole to be deducted, 
g = the distance between gage lines, 

Ss = the pitch or stagger, and 

d = the diameter of the hole. 


Other formulas have been suggested as approximations to the Smith- 
Cochrane formula. One due to C. R. Young® is 


x = 1.50 - 1.25 s/g. (2) 


Cochrane proposed? the formula 


3. “Some Tests Bearing on the Design of Tension Members”, by Edward 
Godfrey, Engineering News, vol. 55, 1906, pp. 488-489. 

4. “Calculating Net Section of Riveted Tension Members, and Fixing Rivet 
Stagger”, by Victor H. Cochrane, Engineering News, voi 59, 1908, p. 465. 

5. “Diagram for Net Section of Riveted Tension Members”, Engineering 
News, vol. 73, 1915, p. 893. 

6. “True Net Sections of Riveted Tension Members”, by C. R. Young, 
Engineering Research Bulletin No. 2, University of Toronto, 1921, 
pp. 233-244, 

7. “Rules for Rivet-Hole Deductions in Tension Members”, by Victor 
H. Cochrane, Engineering News-Record, vol. 89, 1922, pp. 847-848. 
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r=de- (3) 


where r is the width of strip to be deducted. This is the formula in general 
current use, specified for calculating net area for sections with staggered 
holes by the AREA, the AISC, and the AASHO. The following quote is typical: 

“The net width for any chain of holes extending progressively across the 
part shall be obtained by deducting from the gross width the sum of the di- 
ameters of all the holes in the chain and adding, for each gage space in the 
chain, the quantity: 


where s equals the pitch of any two successive holes in the chain, g equals 
the gage of the same holes. 

The net section of the part is obtained from the chain which gives the 
least net width”. 

Tests to failure of drilled steel plates with two and three staggered holes 
were performed by C. R. Young and W. B. Dunbar® who concluded that the 
Smith-Cochrane formula was somewhat too severe. 

Clarence Holleman? reports tests to failure on several aluminum alloys, 
as well as on SAE X4130 steel. Holleman recommends the use of a tension 
effeciency of 90% based on net area for tension members with holes when 
using aluminum alloys. C. H. Stevenson*” reports tests to failure on sheet 
specimens of several aluminum alloys with one and two holes drilled in 
various configurations. Stevenson reports increasing efficiencies as stagger 
is increased, but there are wide variations of efficiency of as much as 13% 
between the alloys tested. In the above papers, efficiency of a tension mem- 
ber is defined as the ratio of ultimate strength based on net area to the ulti- 
mate strength of a solid tensile specimen, expressed in per cent. 


8. “Net Section of Drilled Plates”, by C. R. Young and W. B. Dunbar, 
Engineering Research Pulletin No. 6, University of Toronto, 1926, 
pp. 51-65. 

. “Tension Joints in Aircraft Structures”, by Clarence H. Holleman, 
Journal of the Aeronautical Sciences, vol. 10, 1943, pp. 295-302, 
310. 

“The Effects of Open Holes on the Tensile Strength of some 

Aluminum Alloys”, by C. H. Stevenson, Symposium on Notch Effects 

in High-Strength Aluminum Alloys in the Journal of the Aeronautical 
ciences, vol. 13, 1946, pp. 395-397. 
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H. N. Hill and R. S. Barker!1,12 report tests on specimens cut from 
aluminum and magnesium alloy sheet, the tests being carried to failure. 

The specimens were taken at right angles to the direction of rolling. They 
reported that the specimens containing two staggered holes showed a reduced 
efficiency, based upon the s*/4g rule, the reduction increasing with increas- 
ing ratio of pitch to gage up to a ratio of 1. They also reported that the 
greatest reduction in ultimate strength was found for the most ductile mater- 
ials. 

Unfortunately, as none of these papers reported limit load data, no direct 
comparison can be made with the present tests. 

The hole patterns discussed in this paper were chosen so as to give the 
sgme least net area for two or more rows of holes as for one row using the 
s“/4g rule; see Fig. 2. The configurations were the one row case (1), the 
two row case (1, 2U), the three row cases (1, 2U, 3U) (2L, 1, 2U) termed 
3S and 3L respectively, the five row case (3L, 2L, 1, 2U, 3U), and all nine 
rows. Thus, on the basis of the s2/4g criterion, the sheets with one or 
more rows should all have the same limit load. 


Limit Analysis Theorems 


The limit theorems offer a method of determining upper and lower bounds 
for the limit load of any perfectly plastic material and therefore avoid the 
almost impossible exact plastic analyses which would otherwise be required. 
This method will be described as applied to perfectly plastic materials, ina 
state of plane stress, which obey the maximum shearing stress criterion of 
Tresca. According to this criterion, plastic deformation occurs only when 
the maximum shearing stress has reached a value k, the yield stress in 
pure shear. Perfect plasticity implies that the maximum shearing stress can 
cannot exceed k. If principal stresses 91,09 are employed as Cartesian 


coordinates (a3 = 0), the states of stress with the maximum shearing stress 
k are represented by the points of the hexagon ABCDEF in Fig. 3. 

The principal axes of the plastic strain rate or increment in an isotropic 
material coincide with the principal axes of stress. The principal plastic 
Strain rates are denoted by €y Ey: and € 3” and are in the directions of o> 


>, and Oe, respectively. If the plastic strain rate is represented by a 


vector with principal components € 1’ Eo, and if strain axes are superposed 


on the stress axes, this “vector” is outwardly normal to the hexagon repre- 
senting the state of stress in plastic flow!3 as shown in Fig. 3. At the 


11. “Effect of Open Circular Holes on Tensile Strength and Elongation of 
Sheet Specimens of some Aluminum Alloys”, by H. N. Hill and 
R. S. Barker, NACA Technical Note 1974, 1949. 

12. “Effect of Open Circular Holes on Tensile Strength and Elongation of 
of Sheet Specimens of a Magnesium Alloy”, by R. S. Barker, NACA 
Technical Note 2716, 1952. 

“Some Implications of Work-Hardening and Ideal Plasticity”, by 

D. C. Drucker, Quarterly of Applied Mathematics, vol. 7, 1950, 

pp. 411-418. 

“A More Fundamental Approach to Stress-Strain Relations” by 

D. C. Drucker, Proceedings First U. S. National Congress of Applied 
Mechanics, June 1951, pp. 487-491. 
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vertices, the direction of the vector representing the plastic strain rate may 
lie anywhere between, and including, the directions of the plastic strain rate 
determined by the two adjacent sides of the hexagon. It should be noted that 


the incompressibility condition €, + €) + €,= 0 determines €, when €; 
and Eo are known. 


When a tension member with a row of holes is loaded, the strains are 
purely elastic at first. As the loading is continued, the maximum shearing 
stress at points on the boundary of the holes becomes equal to k in magni- 
tude and plastic flow may take place. Further increase in load causes the 
plastic regions to grow but they are initially completely contained in regions 
in which the strains are still elastic. The strains in the plastic regions, 
therefore, will still be of the same order of magnitude as the elastic strains 
in the surrounding regions. This state, after Prager, is called contained 
plastic flow. At some load, often after regions of plastic flow have reached 
completely across the member, the plastic region is no longer constrained 
by elastic regions. The load at which this occurs is called the limit load, 
and the member is now in a state of impending plastic flow. In the cases 
considered here, large plastic deformations will then take place at no in- 
crease in load. 

The first basic limit analysis theorem* states that the limit load is the 
largest load for which a statically admissible stress field may be found. A 
statically admissible stress field is here defined as one which 


(a) satisfies equilibrium everywhere, 
(b) satisfies the stress boundary conditions, and 
(c) has nowhere a maximum shearing stress larger in magnitude than k. 


Such a stress field may be discontinuous as will be seen later. 

The second basic theorem states that the limit load is the lowest load for 
which a kinematically admissible velocity field can be found. A kinematically 
admissible velocity field is here defined as one in which strain rates are to 
be considered as plastic only and 


(a) the incompressibility condition is satisfied, and 

(b) the rate at which work is done by the external loads or surface trac- 
tions is equal to, or exceeds, the internal rate of dissipation of energy 
for the whole body under consideration. 


This theorem remains valid for velocity fields with surfaces of discontinuity 
through which the velocity component normal to the surfaces varies continu- 
ously while the tangential component is discontinuous. These surfaces may 
be treated as the limiting cases of a transition layer with continuous but 
rapidly varying, tangential velocity components. 

For plane stress and continuous velocity distributions, the internal rate 
of dissipation of energy per unit volume is expressed by 


¢ + oe 


(4) 


The internal rate of dissipation of energy for the entire body is obtained by 
integrating the expression (4) over the volume of the body. Note that the 
internal rate of dissipation of energy for plane stress fields is the scalar 


*For proofs of these theorems, see references 1 and 2. 
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product of the vector from the origin to one of the stress points on the hexa- 
gon of Fig. 3, and the vector representing the plastic strain rate. It may be 
written as 


= wen tel (5) 


where max i€! is the magnitude of the absolutely largest principal compon- 
ent of the plastic strain rate € v €o, or é ,h415, For a surface of discon- 
tinuity, the rate of internal energy dissipation per unit area at any point in 
the surface is 


(6) 


where v is the change in the tangential velocity component across the sur- 
face of discontinuity at the point. The total rate of internal energy dissipation 
for the whole surface is obtained by integrating (6) over the area of the sur- 
face. 

The rate at which work is done by the external loads or surface tractions 
may be computed once a velocity field has been determined. 


Limit Loads for Patterns of Circular Holes 


Consider now the computation of upper and lower bounds of the limit loads 
for each of the six configurations tested: 


(a) One row configuration of five holes. A stress field consisting of 
uniaxial tension in strips lying longitudinally between the holes is statically 
admissible if the tensile stress does not exceed 2k, Fig. 4a. Note that the 
end boundary conditions on the sheet are simply that the sheet is pulled by 
rigid bars rather than that uniform stress is applied. As the sheet widens 
out at the end, however, the two cases are in fact equivalent. Then 


2 = Sd) 2k, (7) 


where 


T* = the limit load for the sheet, 
t the thickness of the sheet, 
W =the width of the sheet, 

d =the diameter of the holes, 


and 
2k = the stress in the strips between the holes. 


. “Limit Design of Reinforcements of Cut-Outs in Slabs”, by P. G. 
Hodge, Jr. and W. Prager, Brown University Report B11-2 to the 
Office of Naval Research, 1951. 

“The Application of Limit Analysis to Punch Indentation Problems” 
by R. T. Shield and D. C. Drucker, to appear in Journal of Applied 
Mechanics. 
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2kvt, (8 ) 
a lower bound for the ratio T*/T,, is 


TH/T = (W - (9) 
One kinematically admissible velocity field for the one row configuration 


is shown in Fig. 4b, half of the sheet sliding with velocity v along a section 
through the holes 45° out of the plane of the sheet. Here, 


Da = kv (10) 


for the whole discontinuity surface, and the total rate of internal energy dissi- 
pation is kv times the area of the surface. Thus, for very thin sheets, the 

total internal energy dissipation rate is 
= (2 t(w - Sd) kv, (11) 


The rate of external work done on the sheet by the tensile load T is 


(12) 


E, = Tv/ 


Equating (11) and (12) gives T which is here an upper bound for the limit 
load T* so that 


T#/T, < (W Sa)/W. (13) 


The same result is obtained with the symmetric discontinuous velocity pattern 
of Fig. 4c. As the upper and lower bounds coincide, the solution for the limit 
load is exact for the maximum shearing stress criterion within the assump- 

tions of the theory for very thin sheets, and 


T#/T, = (W = 0.75 (14) 


for inches and d = 1/2 inch. 
The restriction to sheets of thickness small compared to the diameter of 

hole is clear from Fig. 4. Out of plane sliding on a 45° angle would take 

place on an appreciably larger area than (W - 5d) t V2 for thick plates. 


(b) Nine row configuration. A statically admissible stress field which 
takes out all the holes similar to that for the one row configuration gives 


Te 2 (W - 9d)2k, (15) 


T#/T, 2 (W - = 9.55, (16) 


A simple kinematically admissible velocity field is obtained by letting one 
portion of the sheet slide on two sets of planes perpendicular to the plane of 
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the sheet and through the nine holes on a 45° plane as shown in Fig. 5a. Here 
again, the internal dissipation, Ej, is k times the area of the surface of 
sliding multiplied by the discontinuity velocity v, or 


E, = (% J - (84 + J2 a) kv (17) 


The rate at which work is done by the external force is 


= (18) 
Equation (17) and (21) gives 
Tz < ktLow - (8 + 2) a) (19) 
or 
TH/T, $1 - (bh f2 +1) = (20) 


This result may be improved by changing the angle of sliding slightly, 
Fig. 5b. Although L’ +4L is increased, v cos a likewise gets larger. 
Anangle 8 = 28°, corresponding to a@ = 36° results in 


la 


0.63 (21) 


(c) Two row configuration. The stress field which was used for the nine 
row configuration is statically admissible for the two row configuration. 
Also the velocity field used for the one row configuration is a kinematically 
admissible velocity field for the two row configuration. Thus 


0.55 $ T#/T, $ 0.75 (22) 


In fact, these are bounds for all configurations with two or more rows of holes. 
To improve the lower bound, consider the stress field which is a super- 

position of the stress fields shown in Fig. 6a and 6b. The shaded regions are 

isotropic stress regions, 04 = 09 A detail of the manner in which the stress 


field is taken through the holes is shown in Fig. 7. Fig. 6b incorporates a 
discontinuous stress field! ® the details of which are sketched in Fig. 8a. The 
maximum shearing stress is not increased when the stress field of Fig. 6b 

is superposed on the stress field of Fig. 6a if the angle 6 = 35° and c is 
made as large as possible as shown in Fig. 8a. For the present configuration, 
for y = 19.9°, which minimizes L", andfor @ = 35°, a lower bound on 
T*/T, is 0.70. Therefore, for the two row configuration, 


0.79 < T#/T, £ 0.75. (23) 


(d) Three row configuration 3-S (1, 2U, 3U). Here also a similar stress 
field to that applied to the two row configuration may be used, Fig. 8b. The 
geometry is such that an angle 6 = 33° for the discontinuous portion of the 
stress field is possible without increasing the maximum shear. However, the 
excess Side tensile force is twice as large as for the two row configuration, 
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hence the lower bound of the limit load is correspondingly reduced. For the 
present configuration, 


T#/T, 2 0.66 (24) 


A kinematically admissible discontinuous velocity field is shown in Fig. 9a, 
the upper portion of the sheet moving up with velocity v, the lower portion 
moving down with velocity v, and the two wedges, ABC and DEF, moving 
laterally with velocity v' = v tana. The details of the velocity field along 
BD are shown in Fig. 9b. The upper bound on T*/T for the simplest case 
of a= 45° is 0.73 andfor a= 36° is 0.72, therefore, 


0.66 ¢ T#/T, ¢ 0.72 (25) 


(e) Three row configuration 3-L (2L, 1, 2U). Fig. 8c shows a statically 
admissible field somewhat more elaborate than for the two row case which 
gives 


T#/T. 2 0.67 (26) 


(f) Five row configuration. A velocity field of the type applied to the 
(1, 2U, 3U) 3-S three row configuration is used here, as shown in Fig. 10. 
The upper bound for o = 45° is 0.71 andfor a = 36° is T*/To = 0.69. A 
stress field more elaborate still than in Fig. 8 determines a lower bound of 
0.58 so that 


0.58 T#/7, 0.69. (27) 
The lower bound is not much of an improvement over the 0.55 of Eq. (16). 
Limit Loads for Patterns of Square Hcles 


It is of some interest to note that if the holes in the one row and nine row 
configurations were half inch squares, the limit loads would be exactly 


T:/T,, = 9.75 (28) 
for the one row configuration, and 

T/T. = 0.55 (29) 
for the nine row configuration. Fig. 11 shows the symmetric 45° velocity 


pattern which gives the same result as the simple one-plane slide of the 
type shown in Fig. 5. 


The s*/4g Rule 


A load carrying plate in a riveted joint which is transmitting some or all 
of its load to cover plates or other members differs essentially from the 
perforated plates just described. The rivets do not act simply as plugs but 
take load out of the plate, Fig. 12. Yielding in the sense of reaching a limit 
load may occur across the first rivet holes or by some combination of first 
and second holes, but generally will not involve additional rows because the 
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loaded rivets interfere. One possible pattern is yielding in zones extending 
between the holes as shown in Fig. 13. The lower portion of the plate is 
assumed to move in the direction of T witha velocity v. Each zone in- 
clined at a@ to the direction of T is then stretched at a speed v sin @ and 
sheared at v cos @. For a very thin plate, using the maximum shearing 
stress criterion, this is equivalent to stretching at a speed 


ie. } = 5 (sina+ 1) (30) 


2 
in analogy t th ion f i 
n analogy to 5 + e expression for maximum 


strain in terms of the normal strain € and shearing strain y . An upper 
bound on the contribution to Tf from the length L of the zone of yielding 


is then obtained from 


v sina (z sina e , a 


Tyv = 2kv'Lt 
where t is the thickness of the plate. Thus 
T, < ktL (1 + sin a) (31) 


If, as shown in Fig. 13, the zone of yielding is taken to extend between the 
horizontal diameter extremes at the holes so that 


L sin a = b, = + (32) 


kt fo? + s* +b) 
ktb ( ‘pi + s“/b* + 1) (33) 


ktb (2 + = 


then 


u 


2ktb (1 + = ...) 


which is in a sense an g°/4b rule for the range of s small compared with 

b and is close to the s“/4g specification. This upper ound can be improved 
a Jittle (i.e. lowered) by choosing a better value of a. It appears that 

s“/4g is a reasonable answer for riveted joints where the rivets take load 
out of the plate considered. 


Mises Criterion of Yielding 


If the Mises criterion is employed instead of the maximum shearing stress 
criterion as discussed so far, the yield curve is the circumscribing ellipse to 
the hexagon of Fig. 3. The shearing stress on the planes of sliding discon- 
tinuity is then larger inthe ratio 2: YV3 or 1.155. All computed upper 
bounds for the perforated plates are therefore increased by a little more 
than 15%. Lower bounds remain lower bounds and considerably more elabor- 
ate stress patterns would be needed to raise them. 
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Experimental Program 


Low carbon, hot-rolled steel sheet was used in the first series of tests 
because of its flat yield, and because it was hoped that the mechanical pro- 
perties would be uniform enough so that the tests would not be complicated 
by local variations of these properties. The variation in yield strength was 
so large that several additional series of tests were run on hot-rolled steel 
sheet. 

A 61S-T6 aluminum alloy was chosen to check the results with a material 
which has a much smaller ultimate elongation than hot-rolled steel, and 
which does not have a truly flat yield, Fig. 14. 

All specimens were cut with the long dimension in the direction of rolling 
from sheet approximately 1/16 inch thick. 

The sheet thickness was measured to + 0.0002 of an inch at the intersec- 
tion points of a one inch grid ruled over the central 10 inch long test section 
of each sheet. The readings for each sheet were averaged, the average being 
taken as the sheet thickness. The sheet widths were machined to 10.000 + 
0.005 inches, and the 0.500 inch diameter holes were drilled by means of a 
carefully laid out jig. Lubricated spherical seats and nuts were incorporated 
in the grip system. Two extensometers, designed for simplicity and rugged- 
ness, were mounted one on each edge of the sheet with a 10 inch gage length. 
As a check on the performance of the grips, a dial gage was used to measure 
the relative displacement of the heads of the testing machine. The testing 
machine used showed a maximum error of plus 0% and minus 0.3%, ina 
recent calibration. The general test set-up is shown in Fig. 15. 

The test pieces were pulled to a total extension, in the 10 inch gage length, 
sufficiently large to ensure flow across the specimen, and were not purposely 
carried to fracture. 

Originally, it had been hoped that uniformity of mechanical properties 
would allow the accurate determination of the yield stress of the material in 
each series from three test pieces with no holes. However, uniformity was 
the exception and it was necessary to cut tension test coupons from the test 
pieces themselves, or to use pieces at each end of each large test piece, for 
the purpose of determining the yield stress of that particular test piece. 
Coupons 9.5 inches wide could be obtained from the test specimens because 
yielding in the 10 inch specimens was confined to the region of the holes. In 
the calculations, the yield stress for each large test piece was taken arbi- 
trarily to be the average obtained from the test coupons from that piece. 
This point will be discussed in more detail later. 

Photographs of the flaked scale patterns of representative hot- roiled steel 
test pieces are shown in Fig. 16. It should be noted that these photographs 
show the pattern after the tests were completed. However, the patterns are 
interesting in that the tendency to yield along the diagonal through the holes 
is clearly shown for configurations with two or more rows of holes. Visual 
observation of the one row configuration during tests established that the 
yield occured initially along a transverse section through the holes, as as- 
sumed for the velocity field in the analysis. One sheet was tested with strain 
gages mounted in the transverse section through the holes to establish the 
limit load very accurately and to demonstrate that large strains occur fora 
small bend in the overall load deformation curve. 

A series of curves for specimens from one large sheet of hot-rolled steel 
is shown in Fig. 17. The ordinate is made dimensionless by dividing the 
limit load by the product of the “average” yield stress and the gross area. 
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Starting with the one row configuration, the limit load is marked by hardly 
more than an indentation in the curve; the curves then flatten out increasingly 
at the limit load as the number of rows of holes is increased. For the nine 
row configuration, there are marked upper and lower yields. Fig. 18 isa 
similar set for the 61S-T6 specimens which exhibit the same general pattern 
except, of course, for the absence of an upper yield. Placing plugs in the 
holes has small effect in any case, except for details of shape. 

Figs. 19a and b demonstrate the very large scatter which was obtained for 
the hot-rolled steel sheets. To avoid confusion, the 1, 3, and 9 row cases are 
separated from the 2 and 5 row. All curves are brought in tangent to the 
average elastic line by slight displacements left or right to minimize distor- 
tion of shape. Each limit load selected is indicated by a small arrowhead, 
because in several cases the choice was quite arbitrary. The symbol P 
indicates that the holes were fitted with plugs of almost no clearance. Des- 
pite the scatter, the conclusion is unmistakably that there is an appreciable 
and significant difference between the limit loads for the different configura- 
tions of holes. 

The questions which should be answered are why the data for any one con- 
figuration of holes covers so wide a range and which of the values should agree 
with the theory of limit analysis. The temptation is to take the average value 
and, if this is done, agreement is obtained in principle with the limit analysis 
theory although values are above those calculated with the maximum shearing 
stress criterion. An average for the purpose of comparison of theory and 
experiment is proper, however, if the deviations are due to haphazard factors 
only, such as experimental error. There would seem to be no possibility in 
these tests of instrumental error of appreciable magnitude compared with the 
observed variations so that averaging is suspect. 

Some oversimplified examples explain the difficulty and indicate that the 
lower values of T*/To, the stress ratio, are the ones which should agree with 
an idealized theory. Suppose first that the yield stress is constant across the 
sheet but varies periodically between 25000 and 30000 along the length of the 
sheet. If there are no holes, a yield stress of 25000 will be found experimen- 
tally because the weakest section will govern. A sheet with one row of trans- 
verse holes will exhibit a variation in limit load corresponding to the range of 
25000 to 30000, or over 15 per cent, depending upon whether the line of holes 
is located in a weak or a strong region. If a large number of tests are run, 
the results will scatter as in Fig. 19 and clearly the answer which should 
agree with theory is the lowest limit load obtained, barring experimental er- 
ror. On the other hand, if the yield stress varies in the transverse direction 
but is constant along the length of the sheet, a specimen without holes will 
test at the average value of 27500 and the average of the one row cases should 
be compared with theory. In general, the variation in yield stress will not fol- 
low either pattern and will tend to be somewhat random in addition. Under 
these conditions in most cases there will be some way for the large test cou- 
pons, without holes, to yield at a value close to 25000 while the sheets with 
holes will usually have stronger regions impeding plastic deformation. The 
lower values of limit load, although not necessarily the lowest, should then 
correspond to the theory. 

A variation in thickness has an analogous effect to a variation in yield 
stress but in the experiments did not exceed 2% while the yield stress deter- 
mined by small test coupons in one extreme case varied by more than 25% 
over the entire sheet. An interesting verification of the importance of local 
variations in yield is provided by a comparison of yield stress obtained from 
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the small test coupons, approximately 3 square inches of sheet, and the large 
test coupons of 20 to 60 times this area of sheet. Several cases were obtained 
of small coupons ranging from 25000 to 29000 while the large coupons in the 
same region of the sheet tested between 25000 and 26000. 

A summary of all experimentally determined stress ratios, T*/T, is given 
in Fig. 20 whose points are subject to considerable latitude in their position. 
A glance at Figs. 17-19 shows that the limit load is not defined precisely in 
most experiments with a work hardening material. All indications are never- 
theless that the answers of maximum shear stress theory are reasonable for 
an ideal homogeneous material. However, as usual, a compromise between 
the Mises and the maximum shear stress criteria will probably be better. 

The dimensionless plot enables comparison of the data for aluminum and steel. 
Results for aluminum are in the low range of the results for steel. As the 
aluminum is quite uniform in its properties, this seems to provide further 
confirmation of the influence of variable yield stress in the steel sheet. 


Conclusion 


Both the limit analysis and the tests show clearly that the configuration of 
holes is important in a sheet under tension and must be taken into account 
for an accurate analysis of open or plugged holes. As might be expected this 
is not done properly by the s*/4g rule for riveted joints. The specification 
rule is shown to be based conceptually on an upper bound computation for 
riveted joints in which the rivets take load out of the sheet being analyzed. 

The value of 0.75 for T*/T, is given by the rule for all configurations of 
holes tested. This is appreciably above many of the test results which range 
from 0.63 to 0.92 depending on the number of rows of holes and the uniformity 
of the test specimen. 
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